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Abstract 

In this work we investigate the presence of defect structures in models described by two real 
scalar fields. The coupling between the two fields is inspired on the equations for a multimode 
laser, and the minimum energy trivial configurations are shown to be structurely dependent on the 
parameters of the models. The trial orbit method is then used and several non-trivial analytical 
solutions corresponding to topological solitons are obtained. 
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The study of topological defects is a well established field, particularly for models de- 
scribed by scalar fields [l[ Q]. The simplest topological defect - the kink - arises in theo- 
ries of scalar fields in two-dimensional space-time |3j]. For usual models with spontaneous 
breaking of global symmetry, such defects interpolate between two minima of the potential. 
Important examples in condensed matter physics are the well-known domain walls, which 
separate regions of different magnetization. These defects are essentially classical objects 
with localized and stable distribution of density energy. In the case of two coupled real 
scalar fields, the equations of motion are very hard to solve due to non-trivial nonlineari- 
ties. 
Refs. 

In 1979, Rajaraman proposed a method to solve the pair of equations of motion which 
usually appear in models described by two real scalar fields |7]: it is named the trial orbit 
method, which relies on the search (in a trial way approach) for an appropriate orbit the 
two fields have to obey in the two-dimensional configuration space. Eventually, when one 
tries the right orbit, we can be able to solve the problem analytically. However, since the 
equations of motion are second order differential equations, the task of finding exact solutions 
is very hard and the trial orbit method is not much efficient. 

Some years before - in 1976 - an interesting work [k]] identified an important class of 
models, showing how to reduce the equations of motion to a system of first order differential 
equations. In 1995, this subject was studied by one of us in ref. jll|, that is, the Rajara- 
man's trial orbit method E] was applied for the first order equations obtained within the 
Bogomol'nyi procedure 10| . The use of the trial orbit method for first order differential 
equations was shown to be very efficient and this new procedure allowed us to make inter- 
esting progress, as it is shown in [l2| and in references therein. More recently the use of 



the trial orbit method for models whose equations of motion can be reduced to first order 



differential equations was systematized in 13fl. Other investigations on similar issues have 

nn 

also been done in [14j-[20[, which use distinct procedures and motivations to study two-field 
and other related models. 

In the case of a model with two fields, the kink-like solutions are orbits in the field 
space. In this work, we will further explore the trial orbit method to investigate models 
described by first order equations. Here, although, we construct a class of models inspired 
in a semiclassical theory of multimode laser and use the trial orbit method to find exact 
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solutions that minimize the energy of the field configurations. The results show that, under 
certain conditions on the parameters of the system, several possible solutions connecting 
distinct minima of the models exist. 

We consider a class of models in (1,1) Minkowski space-time dimensions described by the 
relativistic Lagrange density 

c = ^0x^0! + \d^ 2 d^ 2 - vfa, fa) (l) 

where 0i and 2 are the two real scalar fields, and we use the metric such that x° = xq = t 
stands for the time, while x l = —x\ = x represents the spatial coordinate. The notation 
is usual for relativistic theories, with upper (lower) \i standing for contravariant (covariant) 
coordinates. The metric tensor is a diagonal 2x2 matrix, compactly written as = (1,-1). 
The Euler-Lagrange equation, 

o 9C dC 
^<9(<9 M 0) d<p 

leads to the following equations of motion: 

«9 2 0i <9 2 0i dV 

^-^ + w. = (3a) 

d 2 <t> 2 9 2 2 dV 

^-^ + Wz =0 (3b) 

We are interested in kink- like solutions, which are described by static fields - 4>\ = (f>i(x), 
02 = 02 (x) - so that 

d 2 0i _ dV _ d 2 2 _ dV 

dx 2 (90i ' dx 2 90 2 

In general, these equations are very hard to solve, but this task may be simplified if it is 
possible to replace these second order equations by first order differential equations. In 
order to get first order equations, we suppose that the potential is given in terms of another 
function, W = W((f>i,<f>2), as bellow: 

1 fdW\ 2 1 fdW x 2 



In this case, the Bogomol'nyi method allows to argue that the solutions of the first order 
equations 

#i _ dW _ d0 2 _ dW 

dx (90i ' dx 902 



are also solutions of Eqs. (pEj), as it can be easily verified. 

The potential of the above model has zeroes at the singular points of W((fii, fo), and 
this set of singular points forms the vacua manifold of the field theory under investigation. 
Usually, distinct pairs of minima define distinct topological sectors of the model, and the 
solutions of the first order equations are defect structures with an energy cost given by 
E = \AW\, where 



AW = W(M+oo), M+oo)) - W(M-oo), M-oo)) (7) 

with the points (0i(+oo), 2 (+oo)) and (4>i(— oo), 2 ( — oo)) identifying minima in the vacua 
manifold. Since the energy density of the static fields is given by 

, , 1 fd^V 1 fd(f) 2 \ 2 1 fdW\ 2 1 fdW\ 2 

the energy is always positive, and the solutions which obey the first order equations are the 
minimum energy configurations in each topological sector of the model. 
To be specific, let us now consider the superpotential 

W = + - \\i<\>\ - \\22<Pi ~ \\nfy\ (9) 

This choice represents a class of models described by the two sets of parameters: {fix, 1^2}, 
and {An, A22, A12}, the first being mass parameters while the second specifying interactions 
between the two fields. This potential implies the following first order differential equations: 



dx 



- = (flX- Xxx(j>\ - Al202)01 ( 10a ) 

(fi 2 - A 2 i0? - A 22 02)02 (10b) 



where we have set A21 = A i2 . 

The present model represents in reality a family of models which is refereed to it some 
generality. Moreover, there is another specific motivation to adopt it: the system of Eqs. 
(HOP is connected with the semiclassical theory of the laser and can simulate the competition 



22], pp.126-131). It is said 



between two adjacent modes in a cavity above the threshold 
that the laser is at threshold when the pumping rate from the lower state to the upper 
excited state is just sufficient to overcome the cavity loss. In this way, for the particular 
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case of a two-mode laser, within the approximation that the induced transition rate is well 
below the saturation rate, we have (note the resemblance with Eqs. (ITUh)-(ITUb)): 

p2 



E n fi n E n X nn E n ^ ^ XnmEnE^ 



1,2 



(10c) 



Here E n is the time-dependent slow-varying amplitude associated with the mode n, after 
expanding the electric field in the cavity in terms of a complete set of axial modes. With this 
motivation, the parameters /ii and \i2 represent the overall gain, with the condition /ij > 
being necessary to establish the laser oscillation in the mode i [i = 1,2). Furthermore, 
An and A22 are saturation parameters, and one must have \u > for positive population 
inversion of the mode i. The parameter A12 stands for the nonlinear saturation effect on the 
coupling between the two modes. We also have for the two-mode laser A12 ~ A21 > 0. The 
study of competition among modes considers the analysis of the stability of the stationary 
solutions in a phase space diagram of E\ versus Ef , where numerical solutions for arbitrary 
initial conditions reveal the stable and unstable points. It is found that stability of solutions 
is strongly dependent on the parameters, where one can have laser oscillation in just one of 
the modes or a simultaneous oscillation is both modes. 

Our work considers a similar problem. However, instead of investigate (f)\ and <p\ in a 
phase space diagram, we follow another route and make an analysis in connection with the 
field description, searching for analytical description of the fields 0i and 02- To make the 
work as general as possible, let us start considering the vacua manifold, e.g. searching for 
all the possible minimum energy points of the potential, the critical points of W. Initially 
we can count five points (0i, 2 ) of minima: (0, 0), (±0*, 0),(0, ±02) with <p\ = fii/\u and 
02 = \J A*2? ^22 ~~ see Fig. 1. The case where both <pi and 02 are non vanishing can lead to 4 
more points of minima, a continuum of points or no more points, depending on the relation 
between the parameters. We use the first order equations (flOl) to get 



£il - An0i - Ai202 
/i 2 - A 2 10? - A 2 202 








(11a) 
(lib) 



We then define the matrices 



A 



^A A N 

A2I A22 



(a A N 

A*i ^12 



A«>i) = 



Hi A 



22 



^11 A*i 



^2) 



A21 A*2 



V 



(12) 



/ 



^2^12 


— H1X22 


^12 ~~ 


A11A22 




— ^2^11 


^21 — 


A11A22 



We can analyze better the structure of the solutions expressing the former equations in a 
matricial form A$ 2 = /2. 

• For det(A) 7^ we have a formal solution $ 2 = A _1 /7 and the four minima (±0 1; ±02), 
with 

(13) 
(14) 

y ^21 — ^11^22 
See Fig. 1(a). 

• For det(A) = det(A^) = det(A^)) = we have A i2 = ±^X n X 22 . This means 
coalescence between the ellipses represented by Eqs. filial) and flllb|) and we have an infinity 
of solutions. See Fig. 1(c). 

• For det(A) = 0; det(A^ 1 - ) ), det(A^ 2 ^) 7^ there are no solutions satisfying both Eqs. 
filial) and flllbl) and we have a situation of non-touching ellipses. See Fig. 1(b). 

There are other possibilities, which are also shown in Fig. [TJ In the diagrams depicted 
in Fig.l, we show how the minimum energy points change with the signal of the fractions 
A i i/Ai2, 1^2/ X12, y^i/Aii and /U 2 /A 22 . In the following we analyze solutions connecting pairs 
of minima related to the configurations shown in this figure. 

We first deal with the case involving the two crossing lines of minima, as depicted in 
Fig. 1(f). We use equations (TTTT) to get 

J,2 ^ 12 72 72 ^22 12 f-tr\ 

01 = ^ 02; 01 = -T— 02 ( 15 ) 



An™ 1 A 
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These expressions lead to A12 = ±\/AiiA 22 . Now for 1; 2 7^ we have A 12 /An < and 
A22/A21 < 0. Then we have the following choices: (a) if An < =>- Ai 2 > and A 22 < 0, 



or if A 22 < ==>• A21 > and An < 0. In both cases this implies A12 = VA11A22; (b) if 
An > =>- Ai 2 < and A 22 > or if A 22 > =>- A21 < and An > 0. In this case one 



i = -^l =► 0i = ±f^) 02 (16) 



has A12 = — \/AiiA22- For both (a) and (b) cases we will have 

1/4 

*) 

An V^n/ 
Also W(4>i, 2 ) — and there is no kink-like solutions connecting any points in the lines of 
minimum energy. 

The next study concerns the coalesced ellipses of minima, which is depicted in Fig l.(c). 
In this case we have the trivial (0, 0) solution plus a continuum of minima represented by 



the degenerated ellipse. We have W(0,0) = 0, W(±<j>\, 0) = (1/4) (/i?/ An), W{Q,±<f>l) = 
(l/4)(/^2/A 2 2) and W(±0J, 0) = W(0, ±02*)- This means a null energy for all orbits connect- 
ing the coalesced ellipses. The energy of a kink-like structure connecting a point from the 
ellipse and the origin (0,0) is given by E = \ W{0,0) - W(0i,0 2 )| = (V 4 )Oi Ad)- 

The trial orbits method can be used to find an explicit solution for 4>i(x) and 02 (a?) that 
connects (0, 0) — > (0i, 2 )- We try a solution of the form 

0i = A<fi (17) 

To satisfy the minimum energy points one must have A = 0x/0f- Differentiating Eq. (|T71) 
we obtain 

0; = A^f-y 2 =► £ = ^ (is) 

But, considering equations ( TTTJT) we see that this is equivalent to a proportional relation 
among the two ellipses, in the non-degenerated case. We can obtain the B parameter after 
substituting explicitly the equations of the ellipses in Eq. (ITS"]) . This gives 



5) (19) 

and the structure of the orbit depends strongly on the product A11A22, as shown in Fig. [2j 

We now deal with the case of intersecting ellipses, which is depicted in Fig. Ufa). This case 
is very interesting, and it is better to refer to Fig. [3K which shows the general configuration 



for the minimum energy points, where we defined 0i 2 = \J //1/A12 and 02i = \j ^2/^21- 

To obtain the energies of the solutions connecting minima we first consider Eq. 0. We 
have, by symmetry, W(<f>i,<f>2) = ^(01,-02) = W^(— 0i,0 2 ) — W{— 0i, — 02); thus, there 
are no kink-like structure connecting the intersecting points from the ellipses. Also, we have 
W(<f>t,Q) = W(-0[,O) = |/i?/(4A u )| and W(0,(j)* 2 ) = W(Q,-<f>t) = |//|/(4A 22 )|- We studied 
the following cases: 

i) One connection by means of a straight line between (0,0) — > (0, ±0^), with energy 
E x = \W(0,0\ - W(0, ±02)1 = |m|/(4A 2 2) I- This can be found solving the equations of 



motion to obtain 



M«) = o, ^ M = ± ^( 1+tan 2 h( ^ ) ') (20) 



This solution represents a laser operating only on mode 2, where the laser intensity smoothly 
increases from zero to the maximum operating value. By symmetry one can easily find 
similar solutions that connects (0, 0) — > (±0*, 0) where the laser operates only on mode 1. 



ii) We look for solutions that connect (0,0) — > (fa, fa), with energy E% = \W(Q, 0| — 
W(4>i,4>2)\- We try the orbit 

(21) 



Differentiating the orbit and substituting the first order equations leads to 

/ii - \ n <f>l - Ai202 = 5 (/"2 - A 2 10i - A 22 02) ( 22 ) 

Equating the independent coefficients we obtain B = /ii//i2 and the remaining condition 
can be written as 

An - -A 2 il <Pl = (^22 - A12I <& (23) 



/i 2 / VA*2 

K2\ B A2 



$detA<«(iT) = $detA<*>£f (24) 

V 02 / 02 

We then have the following possibilities: 

• det A^ 1 ) = det = 0. This leads to the already analyzed case of coalesced ellipses. 

• detA^ 1 ), detA^ 2 ) ^ 0. This leads to B = 1 and <^detA^ 2 > = ^detA^, which 
means \i\ = /i 2 and An = A22 = A, respectively, with fa = fa. In this case one can obtain 
an expression for <fri (x) by means of substituting the former equation into the equation of 
motion for the field fa(cf. Eq. [10aj ). One finds 



fa(x) = fa(x) = ±y ■ i(l + tanh(/ix)) (25) 

One can see that for x — > —00, (fa, fa) — > (0,0) and for x — > +00, (fa, fa) — > (±fa,±fa). 
This corresponds to a laser operating is both modes with the same intensity, since we have 
for this solution <pi(x) 2 = (p2(x) 2 , which corresponds to Ei(t) 2 = E 2 (t) 2 = I(t). The 
intensity of the i-th mode increases continuously until achieving the maximum value given 

by Imax = M/(Al2 + A). 

iii) For the orbit fa = A§\ + B to connect (0^,0) — > (0,fa), one must have B = 



\/ /xi/An = fa and A = — (X22/ l^2)yJJhjMij = —fa/fa 2 - Then we have the orbit 

4+feY = l (26) 



Deriving the orbit and using equations (|T0|) leads to the consistency conditions 

An = 2A 12 , ^i = 4 + 2^ (27) 
^22 A*2 
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We can show that conditions (|2"T|) are also compatible with the Eqs. (jl~3l) and ffl~4l) that 
define the intersecting points. Also, the condition for existence of the minimum energy points 
02) leads to another constraint on the parameters. We can see in Fig. [3](a) that crossing 
among the ellipses exist only if (f)\ > 2 i and (f)* 2 > 0i2- This leads after using the consistency 
conditions to 2 < /ii//i2 < 2+//i/ /i 2 - This inequality is satisfied only when /ii//U 2 > 2. When 
the ellipses do not cross one another, we must have, as one possibility, that <p\ > 02i and 
02 ^ 012- This, with the conditions (12?|) lead to the condition ^1/^2 > 2 + fii/fi 2 which is an 
impossibility. So, this type of orbit needs the crossing among the ellipses. The phase space 
diagram for this orbit is shown on Fig. [3^b). There one shows that the points (±0i, ±02) 
are unstable. In this way the orbits connect one of these points, for x — > —00, to one of the 
other minimum energy points (0, ±0i) or (±02, 0), for x — > 00. 

As an example we can choose An = 3, A22 = 1/4, = 3,/i2 = 3/4. This means a ratio 
^2 — 4 > 2 and A12 = An/2 = 3/2. This choice corresponds to an initial condition where 
mode 1 (represented by the 0i field) is well above threshold, whereas mode 2 (correspond- 
ing to the 02 field) has a smaller gain. We have minimum energy points (0*,O) = (1,0), 
(0i, ±0 2 ) = (1/2, ±VQ/2), (0, ±0|) = (0, ±v / 3) and an orbit 0i = -(1/3)0^ + 1 connecting 
theses 6 points (three for 2 > and three for 2 < 0). Substituting the orbit in Eq. (j 10 aft , 
we obtain 

st = -5 M - % + 1) (28) 

Integrating the former equation we obtain two solutions that agree with 02 > 0, namely 
((0f (x), 02 (x)) and (fc(x),<f>t(x)), with: 



0f (x) = ~ ( 1 ± 4/1(1 + tanh(3x/4) j (29) 



and 



i 




1(1 + tanh(3x/4) ) (30) 



For (0 X (x), 02~(x)) we have hm^^^^ = 1/2, lim x __ oo 02 = Vd/2 and \xm x ^ +00 4>t = 1, 
lima.^4.00 02 = and the orbit connects (0i,0 2 ) — > (0!,O) as (1/2, VQ/2) — > (1,0). Here we 
have a final state where only mode 1 oscillates. In this regime we can say that oscillation 



in mode 2 was quenched [22J by the oscillation in mode 1. For (0 X (x), 0^ (x)) we have 



lim a; __ oo 1 = l/2,lim^„ oo 2 f = VQ/2 and \im x ^ +OD <f) 1 = 0, linx^+oo 02 = v^, and 
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the orbit connects (^i,0 2 ) — * (0,02*) as (1/2, \/6/2) — > (0, a/3). Now we have the opposite 
regime, where mode 2 absorbs energy continuously and mode 1 decreases until only mode 2 
remains. 

In conclusion, in this work we have used the trial orbit method introduced by Rajaraman 
[?| to investigate first order differential equations which appear when one uses the Bogo- 
mol'nyi approach to study minima energy kink-like solutions |l0|]. We have studied a family 
of models described by two real scalar fields inspired on the theory of two-mode laser. We 
have determined all the minimum energy points in terms of the parameters which specify 
the model. We have found a rich structure of minima, and several analytical solutions of 
the kink-like type, connecting pairs of minima in the field space. In order to correctly map 
our results for (f>i(x) and 02 (x) to the time-dependent problem of the dynamical competition 
between the two modes Ei(t) and Ez(t), we must interpret x from our mathematical solu- 
tions as the physical time t. This is justifiable after comparing Eqs. (HOal) - (llObl) from our 
classical field theory with Eq. (11 Pel) from the semiclassical laser theory. In this way some 
of the exact solutions here studied where used to study phenomenological situations such as 
laser oscillations between two modes in a multi-mode system. 

The authors would like to thank CAPES, CNPq, MCT-CNPq-FAPESQ and MCT-CNPq- 
CT-Infra for partial support. The authors also thank C. Furtado for discussions and M. M. 
Ferreira Jr. for reading the manuscript. 
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ii l /X l2 ,\i 2 /X l2 >0 
\x l /X n ,]x 2 IX 22 < 

det(A ) * 




(I, /X l2 ,\L 2 /X l2 < 
\l l /X n ,\l 2 fk 22 > 
det(A) = 0;det(A* , )^0 




^, /X u <0,H 2 /X l2 >0 
\x l /X n ,ii 2 /X 22 >0 

det(A) = 0;det(A lf >)*0 



l\ 2 >\i 2 /Xj 2 >0 
(x 2 /X u <(x, /X n >0 

det(A) = 0;det(A (,| '' ) )^0 
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\ O) 




H 2 A 12 >ji 1 /X u 
det(A) it 



I ! 2 = V ,V 2 



det(A) = det(A w) ) = 



<l> 2 










\ (/) 



Mi=0 ^=0 
det(A) = det(A* ) ) = 



FIG. 1: Diagrams showing all the possible minimum energy points. 
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FIG. 2: Phase space for \i<i = A11/A12; Af 2 = A11A22 showing orbits connecting the point 
(0,0) to the coalesced ellipses. Orbit (I) is for A11/A22 > 1, (H) for A11/A22 = 1 and (III) is for 
< A11/A22 < 1- 




FIG. 3: (a) Phase space for //1/A12, ^2/^21 > 0; /Ui/Ah, M2/A22 > showing the 9 minimum energy 
points - 4 of them obtained by the intersection points of the ellipses - and some orbits, (b) Phase 
space for An = 2A12 and A11/A22 = 4 + fJ,i/ fj,2, showing the orbits <f>x = ±</>*(— (p^/^ 2 + !)• 
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